Anharmonic vs. relaxational sound damping in glasses: 
I. Brillouin scattering from densified silica 
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This series discusses the origin of sound damping and dispersion in glasses. In particular, we 
address the relative importance of anharmonicity versus thermally activated relaxation. In this first 
article, Brillouin-scattering measurements of permanently densified silica glass are presented. It is 
found that in this case the results are compatible with a model in which damping and dispersion 
are only produced by the anharmonic coupling of the sound waves with thermally excited modes. 
The thermal relaxation time and the unrelaxed velocity are estimated. 

PACS numbers: 78.40.Pg, 63.50,+x, 78.35,+c 



I. INTRODUCTION 



The understanding of the physical processes leading 
to sound attenuation and dispersion in dielectric glasses 
is still far from perfect. At least four distinct causes 
have been discussed: two-level system s E L 1 3, 0, Q , relax- 
ation H, H, > Rayleigh scattering 0, l9L Il0| and anhar- 
monicity. Anharmonicity produces a direct interaction 
of sound with thermal vibrations. One knows that for 
dielectric crystals it is a dominant contribution to the 
damping of acoustic waves . This has been called ei- 
ther "lattice viscosity" llH or "phonon viscosity" [l3j . 
In Akhiezer's treatment |l4. Il5j . the sound wave modu- 
lates the frequencies of the thermally excited modes and 
thereby their effective temperatures. The return to ther- 
mal equilibrium within this vibrational bath, character- 
ized by the mean thermal lifetime of the dominant vibra- 
tions T t h, leads then to energy dissipation and thereby 
to sound damping. This process can be quite local, oc- 
curing between spatially overlapping but not necessarily 
propagating modes 16]. The time Tth just describes the 
return to equilibrium after a perturbation of these vibra- 
tions by an acoustic wave. Described in these terms, the 
Akhiezer mechanism also ought to be active in glasses, 
even though the thermal modes might then be local or 
diffusive ones, leading to what we shall call "network vis- 
cosity" . One of us had already suggested long ago that at 
the high frequencies of Brillouin scattering experiments, 
this Akhiezer contribution should not be neglected in 
vitreous silica [l^. A simulation based on this mech- 
anism was performed for the case of amorphous silicon 
|18| . It showed that the anharmonicity of thermal vibra- 
tions can induce in glasses sound attenuation larger than 
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in the corresponding crystal owing to the enhancement of 
the Griineisen parameters by internal strains. Other de- 
tailed calculations of the anharmonic damping of sound 
in glasses do not seem available in the literature, except 
for extensive predictions of relaxation times based on the 
fracton model [lfj . 

At temperatures above the quantum regime, owing to 
the considerable disorder of glasses, damping can also 
result from unstable structural features relaxing in the 
strain field of the sound wave 0, |(|. These "defects" 
produce what is called "thermally activated relaxation" . 
The damping depends then on the defect density, on the 
strength of their coupling to strain, and on the distri- 
bution of their relaxation times r. These times are of- 
ten much longer than typical values of r t h- Hence, this 
type of damping appears dominant in many glasses at 
sonic and ultrasonic frequencies. It produces peaks in the 
temperature dependence of the sound attenuation which 
have been well described by phenomenology [(| , although 
the exact nature of the defects remains most often un- 
known. A puzzling aspect is the relative "universality" of 
these relaxation phenomena, as reviewed in |20j ] . Another 
non-trivial aspect is that within the soft-potential model 
the relaxational features at higher temperatures and the 
properties associated with two-level systems, at He tem- 
peratures and below, can be quite succesfully described 
using a single picture 0. 

Another possible damping mechanism is scattering of 
the sound waves by static density or elastic constant fluc- 
tuations. This effect which grows with the fourth power 
of the frequency was described long ago by Lord Rayleigh 
8]. Although this mechanism has been invoked to ex- 
plain the plateau in the thermal conductivity of glasses, 
or the "end" of acoustic branches [lfjl l2ll l22l |23| , it has 
been known for quite some time that it is difficult to con- 
struct a model for glasses that produces strong Rayleigh 
scattering [24| . It has been recently reemphasized that 
the Rayleigh scattering of sound is generally too weak 
for it to play any appreciable role at frequencies below at 
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least 1 THz [2^, [2^. Thus, this will not be of importance 
to the present series of papers as our considerations will 
be restricted to frequencies below ~ 100 GHz. Similarly, 
another very strong damping mechanism is the resonant 
interaction -or hybridization- of the sound with other 
low frequency modes, which can lead to "boson peaks" 
[27l l28l | . In the case of densified silica, this effect has been 
observed to start around 1 THz and to grow then with 
a high power of the sound frequency (23, H3|- Again, 
this will not be of relevance here. Our interest will in 
fact mainly be in the relative strength of the anharmonic 
and relaxation sources of damping. At very low temper- 
atures, below liquid He, the quantum aspect of two-level 
systems becomes crucial 0, but anharmonic damping 
tends to zero owing to the weak population of the ther- 
mal bath. Very low temperatures will not be discussed 
here. 

Compared to the very abundant ultrasonic results 
there exist relatively few damping data at frequencies 
between 1 GHz and 1 THz. The "hypersonic" regime ob- 
served with optical Brillouin scattering falls in the middle 
of this range. The longitudinal acoustic (LA) waves are 
active in the backscattering geometry. To fix ideas, in 
the case of silica and using the green argon-laser line, the 
observed hypersonic frequency is then = D,/2tt ~ 34 
GHz. Explanations that are only based on thermally ac- 
tivated relaxations generally meet difficulties in trying to 
account simultaneously for hypersonic and ultrasonic at- 
tenuation results, e.g. in [32]. However, it was recently 
claimed that for silica, "within a factor of two" , the re- 
sults over the entire frequency range are well accounted 
for by relaxation only |3Jj . A main purpose of the present 
series is to clarify this issue. This becomes all the more 
important that data are now also obtained at frequen- 
cies above ~ 250 GHz, using x-ray Brillouin scattering 
[23, [13 ■ -"- n that context, one should beware of sweeping 
interpolations of the damping behavior over great many 
decades in frequency (34|. As explained in the second 
paper of this series (II) , the damping contribution of an- 
harmonicity in the hypersonic regime increases faster in f2 
than that of relaxations. At sufficiently high frequencies 
and temperatures T, this can lead, over certain frequency 
and temperature ranges, to a dominance of anharmonic- 
ity. Indeed, the thermal phonon contribution increases 
with the thermal population and thus with some power 
of T, whereas the effect of the relaxing defects usually 
saturates. 

In order to estimate the anharmonic contribution to 
sound dispersion and damping, it is of interest to con- 
sider several varieties of silica glasses besides usual vitre- 
ous silica, u-Si02- For example, it is already known that 
irradiation of u-Si02 by fast neutrons increases slightly 
the glass density p, from 2.20 to ~ 2.26 g/cm 3 |35| . 
but thereby increases appreciably the hypersonic veloc- 
ity while it reduces substantially the hypersonic attenu- 
ation (3|J. Similarly, it has already been shown that in 
permanently densified silica glass both the relaxation at 
sonic frequencies [37J and the hypersonic linewidth [3^ 



are very strongly reduced. The latter experiments were 
performed on a sample of mean density p ~ 2.45 g/cm 3 
which unfortunately was quite inhomogeneous. It is thus 
of considerable interest to obtain detailed T-dependent 
Brillouin data on a much better sample of well densi- 
fied silica, <i-Si02, with p ~ 2.60 g/cm 3 . The results of 
such an experiment are presented and discussed in this 
first paper. They confirm that the hypersonic damping 
in <i-Si02 is indeed entirely dominated by anharmonic- 
ity. This is of particular interest as it reveals for the first 
time what the shape and strength of a pure anharmonic 
contribution can be in a glass. 

Two quantities will be of main interest in these papers, 
the damping of sound and its velocity dispersion. Damp- 
ing can be expressed in several ways. It is manifested 
in a Brillouin scattering measurement by a linewidth 
T, where V /2tt is defined as the deconvoluted frequency 
half-width at half-maximum of a damped harmonic os- 
cillator response adjusted to the Brillouin peak. On the 
other hand, in ultrasonics the damping is generally mea- 
sured as an energy mean-free path I, or as its inverse 
l~ l = a = 2T/v, where v is the sound velocity and T 
is in rad/s. Traditionally, the attenuation constant a 
is often expressed in dB/cm, in which case its numer- 
ical value must be multiplied by 10 In 10 to convert it 
to m . To compare these various results over a broad 
range of frequencies, we prefer to use the internal friction 
Q^ 1 which is the inverse of the quality factor Q = Q/2T. 
Thus, Q^ 1 = £~ 1 v/fl. This is particularly convenient 
for graphical presentation as near a relaxation peak Q^ 1 
depends little on O, this over many decades of il. The 
dispersion in v(il, T) is usually quite weak, and thus it 
is indicated to present a differential quantity, Sv/v. In 
defining this quantity a suitable reference value vq must 
be selected, as Sv/v really means (v(Q,T) — v )/vq. A 
rather common choice is to take for the value of v(il, T) 
at the lowest measured temperature, usually near liquid- 
He temperatures. One should be aware that an arbitrary 
choice of vq modifies Sv/v practically by an additive con- 
stant, since Sv <C v. Further, if Q _1 is small compared 
to 1, which is always true here, one easily shows that 
Q" 1 and —2Sv/v are Kramers-Kronig transforms of each 
other, 

-2Sv(n, T)/v = - P / +0 ° 9 ~ 1(3: ' T) dx , (1) 

7T J_ OQ X-il 

with the related relation for the inverse transform. As 
explained in |39|. this strictly requires that Sv = v — 
Voo, where Voo = v(fl — > 00, T) is the high-frequency 
limit of v, which can be called the unrelaxed velocity. 
It is conceivable that "static" structural changes might 
occur in function of T, and that these could produce a 
temperature dependent Voo(T). 

As discussed in if Q _1 is given by a single Debye 
relaxation, 

Q _1 (fi,T) = Afir/(1 + n 2 r 2 ) , (2a) 
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where A and r are independent of frequency but might 
depend on T, it follows from (1) that 

-2Sv(Q,T)/v = A/{1 + Q 2 t 2 ) . (2b) 

Finally, concerning the various independent contribu- 
tions to r, and thus to we shall adopt the equiv- 
alent of "Matthiessen's rule" [41j that these are simply 
additive. Thus, if Q~, is produced by relaxation and 
Q&nh by anharmonicity, the total Q _1 will simply be 
Qrci + Qanh- ^ should also be noted that the right- 
hand side of (1) only describes velocity changes associ- 
ated with frequency dependent damping, while "static" 
velocity changes (t>oo(T) — vo)/vq might result from In- 
dependent structural effects. 

The paper is organized as follows. In Section II, the 
experimental aspects are described, including the spec- 
trometer and the characterization of the sample. In Sec- 
tion III, the new Brillouin scattering data obtained on a 
well densified sample of d-Si02 are presented. These are 
then analyzed in Section IV. A discussion concludes the 
paper in Section V. 

II. EXPERIMENTAL ASPECTS 

The experiments reported here are particularly de- 
manding both in terms of sample quality and spectrome- 
ter resolution. The Brillouin frequency shift vb produced 
by light scattering from an acoustic wave of velocity v in 
an isotropic medium is given by 

/Tin 

vb = -F sin (6/2) , (3) 

where n is the refractive index of the sample, Al the 
vacuum wavelength of the exciting laser, and 9 the scat- 
tering angle. The scattered light must be collected over 
a finite aperture A9. Hence, if one is interested in small 
linewidths - which is eminently the case here - the broad- 
ening produced by A9, proportional to cos((9/2) AO, must 
be minimized by working with the smallest acceptable A9 
and near backscattering, i.e. for 9 ~ 180°. In isotropic 
media the transverse acoustic modes are not active in 
backscattering. This is why the measurements are re- 
stricted here to the LA-mode. 

Densification has an enormous effect on v which at 
room T changes from ~ 5950 m/s for v-SiC>2 to ~ 7050 
m/s for d-Si02 with p = 2.63 g/cm 3 . Thus, near 
backscattering and at Al = 5145 A, vb changes from 
~ 34 GHz in v-Si0 2 to ~ 42 GHz in d-Si0 2 . On the 
other hand, densification also reduces considerably the 
sound attenuation and hence the half-width of the Bril- 
louin line. In v-Si02 at room T a half-width A^b — 75 
MHz is observed [jjj, whereas it falls below 20 MHz in 
c?-SiC>2. Note that Ai>b would equal T/2ir - the quantity 
of main interest here - in absence of parasitic sources of 
broadening. In fact, these narrow linewidths can only be 
measured on samples that are very homogeneous in their 



density. It is easily estimated from the above figures that 
an inhomogeneity within the scattering volume of only 
1% in density would lead to a broadening of the Brillouin 
signal by about 400 MHz. This value would be totally 
unacceptable when the instrinsic width of interest is only 
a few MHz. Furthermore, such small linewidths are usu- 
ally not measurable in commercial Brillouin instruments. 
The Brillouin spectrometer used in these measurements 
is briefly described in Subsection II A. 

Permanently densified silica, d-Si02, is obtained by 
submitting v-Si02 to high pressure at elevated T |42|. 
Our scries of samples was prepared from short cylin- 
ders of silica that were submitted to ~ 8 GPa of quasi- 
hydrostatic pressure at T ~ 700 °C. Different densities 
were achieved by varying the duration of the treatment. 
We used samples with initial p values of 2.31, 2.46, and 
2.63 g/cm 3 , besides w-Si02 with p = 2.20 g/cm 3 . A typ- 
ical raw sample had a diameter of 4 mm and a length of 
8 mm. Although the samples are often cracked, those at 
the highest densities can be clear and transparent. Their 
surfaces are smooth but not flat. Thus, to perform Bril- 
louin scattering at room T, the samples can be placed 
in an index matching fluid. This is not feasible at either 
cryogenic or elevated temperatures. Hence, a selected 
sample of the highest density had to be cut and polished 
with flat faces. It was then carefully characterized to 
identify a suitable region to perform these experiments. 
This is described in Subsection II B. 



A. High resolution spectrometer 

Three main criteria are to be satisfied by a spectrom- 
eter suitable for these Brillouin measurements: contrast, 
resolution, and accuracy. A very high contrast between 
the Stokes shifted frequency Ug = Vi, — vb and the excit- 
ing laser frequency z/l can be achieved by multiple passes 
through a planar Fabry Perot (PFP) interferometer |43j . 
However, the resolution of such a PFP is insufficient to 
measure Ai^b- It is thus placed in series with a spherical 
Fabry Perot (SFP) which is the scanned instrument. The 
high accuracy on vb is achieved by producing a frequency 
adjustable reference signal vm — with an electro-optic 
modulator working near vb- In practice, this reference is 
also used to dynamically adjust the PFP which works as 
a filter only transmiting the Stokes Brillouin line. 

We used a modern version of the instrument described 
in 44]. The PFP has a spacing of 2 mm and it is used 
in four passes. The SFP has a free spectral range of 1.5 
GHz and a finesse of ~ 40. A typical spectrum obtained 
in scanning just the SFP is illustrated in Fig. 1. One 
recognizes two orders of the elastic line at z/l , the Stokes 
Brillouin signal at z/g, and the modulation signal at pm- 
All these lines are at a rather high interference order. A 
very high accuracy on us is obtained as only the fractional 
order separating it from um needs to be determined, while 
z/m itself is known to the extremely high accuracy of the 
frequency generator. This provides for a determination 
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FIG. 2: Hypersonic velocity measured in d-SiCh in function 
of the sample density. The line is a guide to the eye showing 
the initial slope. 



FIG. 1: Signal obtained in scanning a full order of the spher- 
ical Fabry- Perot. One recognizes two elastic peaks marked 
ul, the Stokes shifted peak us, and the modulation signal 
um- The latter is stopped during the acquisition of the Stokes 
signal. 

of ub to better than 2 MHz absolute accuracy provided 
the Brillouin signal is appropriately intense and symmet- 
ric. The instrumental half-width of this SFP is ~ 20 
MHz. Using photocounting and numerical deconvolution 
the half- width of a sufficiently intense damped harmonic 
oscillator (DHO) Brillouin line can then be determined 
with ±2 to 3 MHz precision. 

B. Sample characterization 

The mean densities of the various samples were care- 
fully determined by the flotation method. Their refrac- 
tive index at the sodium D lines (A = 5893 A) was mea- 
sured with an Abbe refractometer. This gave results in 
good agreement with a previous determination |45| . The 
corresponding indices for Al = 5145 A were extrapolated. 
Their values vary linearly with p and are given by 

n(5145A) = 1.022 + 0.200 p (g/cm 3 ) , (4) 

to 10 -3 accuracy. 

Figure 2 shows the hypersonic velocities measured at 
room temperature on the four original densities, plus one 
at an additional density of 2.60 g/cm 3 (see below). The 
samples of intermediate densities (2.31 and 2.46) were 
much cracked internally, which produced a spread of ve- 
locities that translated into broader and distorted Bril- 
louin lines. This did not prevent a sufficiently accurate 
determination of a mean velocity, with a ±10 m/s error 
bar which is still small on the scale of Fig. 2. The impor- 
tant point is that v is a strong function of p. One finds 



dv/dp = 1.87 x 10 3 in m/s per g/cm 3 , at small densifi- 
cations. As seen in Fig. 2, this slope steepens towards 
higher densities. 

Brillouin scattering is a rather local probe. In the 
backward geometry used here, different small areas of 
a sample cross-section can be probed. The measurement 
of ub is then a test for the sample homogeneity. This 
revealed that one of the samples of highest density was 
rather homogeneous, with ub ranging from 42.32 to 42.35 
GHz across its section. This piece was then selected for 
cutting and polishing, to obtain flat faces in order to 
be able to perform the Brillouin measurements in func- 
tion of T . Unfortunately, this mechanical process intro- 
duced an appreciable distribution of frequencies, to such 
an extent that the Brillouin lines became quite asymmet- 
ric. This sample was then used for Brillouin scattering 
with x-rays in |2^| , and subsequently in and |2^| . It 
should be remarked that these inhomogeneities do not 
appreciably affect the x-ray Brillouin scattering results. 
The latter experiments are indeed performed with four 
orders of magnitude less frequency resolution than here, 
the error bars on r/27r being then at best ±30 GHz [29l| . 
During the x-ray measurements the sample was main- 
tained at an elevated temperature for a long time in order 
to increase the signal by the thermal population of the 
acoustic mode. We used T = 575 K [23], a temperature 
which is too low to appreciably relax the densification 
over the long duration of the x-ray experiment. Check- 
ing the optical Brillouin spectrum after the first x-ray 
experiment we discovered that the sample had become 
locally much more homogeneous, presumably owing to 
the thermal treatment. Nice symmetric and narrow Bril- 
louin lines were then observed. At room T, one end of 
the sample had a velocity indicating p ~ 2.60 g/cm 3 , 
while near the opposite end we found p ~ 2.57. The 
corresponding difference in ub is of the order of 1 GHz. 



The former region had a linewidth Aj^b = 25 MHz, while 
the latter had Av^ — 40 MHz. These two regions were 
separated by a less homogeneous one, with asymmetric 
Brillouin lines and a higher width. For the measurements 
presented in Section III, the region of density 2.60 and 
with the narrower line was selected. It should be remem- 
bered that part of the observed broadening can easily 
arise from small remaining inhomogeneities. 



III. BRILLOUIN SCATTERING RESULTS 

Figure 3 presents the Brillouin results obtained on the 
sample region with p = 2.60 as explained in Subsection II 
B. For the low-T measurements, the sample was placed 
in a He- flow optical cryostat. The temperature of the gas 
flow is stabilized to much better than 1 K, and the sample 
temperature is read with a silicon diode. For the high- 
T measurements, the sample was mounted in an optical 
oven stabilized to ±1 K. The value of 9, the scattering 
angle internal to the sample, was taken at 178.1° for these 
measurements in order to minimize the collection angle 
broadening. The measured spectra were fitted to a DHO 
response function of frequency i>b and half- width A;/b, 
convoluted with the experimental instrumental function. 

The values v& are shown in Fig. 3a. The error bars 
on these values are of the order of the size of the points. 
The values Aub are shown in Fig. 3b with error bars 
that are given by the Marquardt-Levenberg least-square 
fitting algorithm. Care was taken to keep the laser beam 
on the same sample spot during the several days needed 
for measurements. It was however not possible to recover 
exactly the same position on the sample after its transfer 
from the cryostat to the oven. If residual inhomogeneous 
broadening is present - as we suspect - its amount is 
likely to vary from place to place. This could explain the 
difference - beyond the size of the error bars - between 
the level of A^b observed in the cryostat and that in the 
oven. 

It is now of interest to briefly compare the results ob- 
tained on d-SiC"2 to these in i;-SiC>2. Fig. 4 illustrates 
data obtained long ago on i>-Si02 [47J. The velocities 
below room-T are compared in Fig. 4a. The values of v 
for d-SiC>2 are derived from the data shown in Fig. 3a 
using Eq. (3). As T is raised from liq. He tempera- 
tures, one clearly observes a dip of v in v-SiC^- Within 
the accuracy of our experiment, such a dip is totally ab- 
sent in d-SiC-2. It is well known that the dip is produced 
by the T-dependence of thermally activated relaxational 
motions. This will further be discussed in paper (II). 
Its absence in d-Si02 demonstrates that thermally acti- 
vated relaxation becomes negligible after densification. A 
strong decrease of these relaxations compared to 1>-Si02 
has already been reported for the less densified sample of 
rf-Si02 that was examined with both Brillouin scattering 
|38j and vibrating reed |37| measurements. In that sam- 
ple there remained a clear dip in v, but the relaxation 
peak observed with vibrating reed had decreased by a 
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FIG. 3: The measured Brillouin frequency shift vb (a), and 
the half- width Ai/B (b), in function of the temperature T for 
the sample of density p = 2.60 g/cm 3 . Different symbols are 
used for data taken in the cryostat and in the oven. Although 
the raw data in (b) might suggest a dip around room tem- 
perature, it is more likely that this effect is related to sample 
inhomogeneities, as explained in the text. 



factor of about six compared to u-Si02- We conclude 
that in our case, the dip in v being absent, the relaxation 
contribution to the damping must be totally negligible. 



Fig. 4b shows the corresponding half-widths. A loga- 
rithmic scale is used here to emphasize the low-T region 
and also because the damping in w-Si02 is quite a bit 
larger than in e?-SiC>2. As discussed in paper (II), the v- 
S1O2 values result from a combination of relaxation and 
anharmonicity. Below 10 K, they fall however below the 
A^b values observed in d-Si02, which would be surpris- 
ing were it not for the inhomogeneity of the c?-Si02 sam- 
ple dicussed above. Hence, we attribute the level of 8 
MHz, illustrated by the dashed line, to the residual inho- 
mogeneous background. The corresponding level for the 
data measured in the oven is estimated to be about half 
of that, i.e. ~ 4 MHz. 
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FIG. 4: Comparison of d-Si02 and u-Si02 showing the ve- 
locities (a) and the linewidths (b) below room temperature. 
The v-Si02 data was acquired long ago [17] with a similar 
spectrometer but at a laser wavelength of 4880 A and with- 
out electro-optic modulated reference. The solid lines result 
from calculations explained in the text. The dashed line in 
(b) indicates the estimated inhomogeneous background of the 
d-Si02 sample. 



IV. NETWORK VISCOSITY 

The data shown in Fig. 3b, corrected for the estimated 
inhomogeneous background, is used to calculate the in- 
ternal friction Q^ 1 = 2T/fi shown in Fig. 5. As ex- 
plained above, this result should be interpreted in terms 
of the Q~n h produced by anharmonicity. The theory 
of lattice viscosity was developed for crystals, in which 
case the calculation of correlation functions |4g or the 
sums over modes ^2 EZJ could be carried out quite far. 
For glasses, only the fracton model has been pursued to 
the point of making analytic predictions for the acoustic 
linewidth 0,|5(]]. However, the development described 
in , in particular the one leading to Eq. (131) of that 
review, is essentially based on the standard equations of 
viscoelasticity supplemented by a simplifying assumption 
regarding the details of the collision term in the Boltz- 
mann equation. A fully parallel reasoning can be made 
for glasses. Hence, for fir t h <C 1, similarly to the phonon 
viscosity term of |l3l . we posit that network viscosity in 
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FIG. 5: The internal friction of the d-Si02 sample. The cal- 
culated line is explained in the text. 



glasses leads to the internal friction 
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(5a) 



(56) 



Here C v is the specific heat per unit volume, v is the 
velocity of the sound wave, vn is the Debye velocity, and 
7 2 is a mean squared average Gruneisen parameter for 
the thermally excited modes 0, 0. We neglect the 
thermal conduction term in Eq. (131) of Experience 
shows that the latter is already considerably smaller than 
the phonon- viscosity one in crystals [5lj . This should 
be even more so in glasses owing to the much reduced 
thermal conduction |18| . The strongest assumption in (5) 
is presumably the use of a single relaxation time Tth(T). 
A second assumption is to take 7 2 constant. 

As the temperature is lowered, one expects that r t h 
increases, so that the condition firth <C 1 eventually be- 
comes violated. We use the Ansatz that Q~ n \ is then 
given by (2a), where r = rth and A is given by (5b), 



Qanh=^W(l+fi 2 r t 2 h ), 



(6) 



As pointed out in [48j], this amounts to introducing an 
exponential decay in the time-dependent relaxation func- 
tion. This is preferable to the expression (4.15) proposed 
in ^ij, as the latter does not have a finite Kramers- 
Kronig transform. 

It is now of interest to estimate rth from the data shown 
in Fig. 5. This is done using (6), in which all quantities 
are quite well known except for the value of j 2 that en- 
ters A according to (56) . The other non-trivial quantities 
in A are: 1) the specific heat C V (T) for which we use the 
data from [5^. complemented at high-T with the spe- 
cific heat of i>-Si02 appropriately scaled to account for 



7 




Temperature (K) 



FIG. 6: The thermal relaxation time obtained from the data 
in Fig. 5 using eq. (5) and (6) as explained in the text. The 
dashed line shows the 1/T behavior at high T. The solid line 
is a guide to the eye explained in the text. 



the higher density of d-Si02; 2) the Debye velocity for 
which we use Vfy — 0.322 v 3 , a numerical coefficient quite 
universal for glasses and taken from v-SiC^. With the 
above, one calculates A/7 2 , which is a quantity that in- 
creases rapidly with T at low T. One divides then Q _1 
by that quantity, to obtain a result which according to 
(6) is equal to 7 2 fir t h/(l + fi 2 r 2 h ). One notes that in 
function of T the latter is maximum at firth = U an d 
that the maximum value is simply 7 2 /2. Unfortunately 
the noise in Q _1 /(A/7 2 ) becomes much too high at low 
temperatures owing to the smallness of A/7 2 , so that a 
clear maximum cannot be identified. One notes however 
that the maximum in Q~u/A is not a maximum in 
as A is a rapidly growing function of T. At firth = 1, 
there is rather a step-like increase in <5 _1 (T). Such a 
feature is clearly recognizable between 25 and 35 K in 
the semilog presentation of Fig. 5. This suggests that 
fir t h reaches 1 around 30 K. An extrapolation of Q~^ h /A 
from higher T indicates then that 7 2 w 8 is an appropri- 
ate value, with an uncertainty smaller than 20 %. The 
resulting value 7 ~ 2.8 is reasonable, especially owing to 
the observations made in 0, [5j| . To complete the task, 
one can then solve for fir t h the quadratic equation (6). 
For each data point, the lowest root must be kept for 
firth < 1 and the highest one for firth > 1- The noise 
in our data produces some complex roots for T < 40 K. 
These have been rejected. Dividing the roots by fi, the 
values of rth shown in Fig. 6 are then obtained. The 
uncertainty in 7 2 produces a similar uncertainty in the 
absolute scale of rth- This is of no consequence for the 
following discussion. 

We observe that r t h increases rapidly on cooling. It 
varies approximately in 1/T from ~ 150 to 600 K, and 
progressively faster below ~ 100 K. This is illustrated by 



the dashed line of slope -1 in Fig. 6. To produce the 
guides to the eye used in Figs. 4 to 6, we first adjusted 
i~th{T) to a sum of inverse powers of T, obtaining the 
solid line shown in Fig. 6. We also adjusted v(T) to 
i>o[l — a/(e e / T — 1)] which is derived from Wachtman's 
equation as discussed in 0. This gives vo — 6853 m/s, 
a = -1.052 x 10~ 2 , and O = 328 K, corresponding to the 
line shown in Fig. 4a which provides an excellent fit of v 
up to the highest measured data point at T near 600 K. 
Having v(T), we can calculate fi(T), and then Q~^ h (T) 
and Avq(T) as shown in Figs. 5 and 4b. The sharpness 
of the step near 30 K strongly depends on the assumed 
T-dependence of rth(T) in that region. Its exact position 
also depends on the value selected for 7 2 . At this stage, 
the limited quality of the data would not justify attempt- 
ing any further improvement on these adjustments. The 
main difficulty obviously arises from the inhomogeneous 
background contribution to Ai^b- 

We can now use Eq. (2) to calculate the velocity 
change (<5u/i>)anh associated with Q~^ h . As indicated 
in the Introduction, the result of this calculation is 
then referred to the unrelaxed velocity Woo(T) whose In- 
dependence can have a structural origin. Thus one has 

v-VcoiT) A/2 

Voo (T) ' l + fi 2 T 2 h ' ^ 1 

where v is the Brillouin value measured at frequency fi. 
All velocity changes being relatively small, Voo(T) can be 
replaced by «oo(0) = Urj i n the denominator on the left- 
hand side. Further, A being anyway very small when 
firth > 1, firth can be neglected in the denominator on 
the right-hand side of (7), so that finally 

Voc{T) -v v - v A 

— r 77 ■ (° ) 

v v 2 

We remark that the size of the "correction" term A/2 is 
directly proportional to the value selected for 7 2 . Figure 
7 shows (v — vo)/vo and (voc — vo)/vq. Both exhibit an 
anomalous T-dependence, increasing as T increases, but 
this anomaly becomes quite a bit stronger for the unre- 
laxed velocity (woo — vo)/vq. 

V. DISCUSSION 

The above results suggest that in d-Si02, contrary to 
all other investigated glasses, the damping of sound ob- 
served in Brillouin scattering is purely produced by the 
network viscosity, i.e. by the anharmonic interaction of 
collective modes rather than by the relaxation of local 
"defects". Two features deserve additional discussion: 
1) the temperature dependence of the relaxation time of 
the dominant thermal modes, r t h(T); 2) that of the un- 
relaxed sound velocity Voo(T). 

The particular behavior of r t h (T) , which decays ~ 1 /T 
above ~ 150 K, probably reflects the fact that in a glass, 
as opposed to crystals, the thermal modes are mostly 
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velocity changes. The line is the adjustment with Wachtman's 
equation (54|. 



diffusive or non-propagating. By comparison, in crys- 
tal quartz for example there is a large range of T where 
Tth decreases with 1/T 2 (l6||. To understand the nature 
of this difference, we consider the only available calcu- 
lation of anharmonic lifetimes in disordered systems. It 
is based on the fracton model [l^. This model accounts 
for the non-propagating nature of the modes. As pointed 
out by the authors of fractons have been used to 
"obtain explicit results" , but "the approach would ap- 
ply to any localized vibrational states" . We believe this 
remark is in particular true for what concerns the func- 
tional dependence in T of the relaxation times. Indeed, 
the dependence found in [l9| does not involve any of the 
various fractal dimensions, but just the crossover energy 
hcu co above which the modes become non-propagating. 
The T-dependence results from the thermal population 
of the modes and their spatial overlap. It is presum- 
ably of quite general applicability to such situations, in- 
dependently from the fractal model. The crossover en- 
ergy is well known for d-Si02 as it has been carefully 
determined by x-ray inelastic scattering measurements 
[23L I29L |4^ | . It is located at about 9 meV, or correspond- 
ingly at T co = huj co /kB — 110 K. The prediction of [lj| is 
that all relaxation times vary like r oc 1/T for T > T co , 
whereas they diverge with r oc exp(T co /T) at lower tem- 
peratures T < T co . We note that, qualitatively, this is in 
perfect agreement with our observations. 

Further, the thermal conductivity k in the Debye 
model can be calculated from tq, using the standard ki- 
netic equation k = C v WQTth/3 |41| . Doing this, we find a 
calculated n which is considerably larger than the mea- 
sured one [55| ■ Its T-dependence also disagrees with the 
observations, as the calculation gives a k(T) that de- 
creases with increasing T at high T, like in a crystal, 
contrary to the observed increase |55l beyond the typi- 
cal thermal conductivity plateau [56( common to glasses. 



The reason for this is of course that the modes contribut- 
ing to C v are at best diffusive so that the kinetic equation 
is not applicable. This provides a second and indepen- 
dent indication that the thermal modes are essentially 
non-propagating. 

Regarding the sound velocity, we first remark that the 
network viscosity has the usual effect of decreasing the 
velocity from Voo to v as T increases. This is evident 
from a comparison of the two curves illustrated in Fig. 7. 
What is abnormal is that the unrelaxed velocity i>oo(T) 
increases with T. This anomaly is common to glasses 
that contain tetrahedral building blocks, as shown ex- 
perimentally for SiC>2, Ge02, BeF2, and zinc phosphate 
in [57]. An earlier proposal that this might be a man- 
ifestation of large structural inhomogeneities special to 
refractory glasses [H^ does not seem supported by the 
fact that the effect is particularly strong in BeF2 which 
has a low glass transition temperature T g . A more likely 
explanation might be that there is a specific structural 
feature in SiCVlike glasses which leads to a T-dependent 
reversible structural change as shown in recent simula- 
tions It was found in these that the angular posi- 
tion of the Si-O-Si plane around the Si-Si direction has 
distinct equilibrium values separated by 90° jumps. The 
latter would produce in the glass a gradual transition sim- 
ilar in nature to the a — (3 transformation in cristobalite. 
One might ask why there is no corresponding contribu- 
tion to Q _1 owing to such a transition-like mechanism. 
A possible explanation it that the transition being dis- 
tributed in T, most of the volume is always very far from 
ilr ~ 1 so that the additional damping might be much 
too weak to be observable. 

Another question of interest is whether the unrelaxed 
velocity Voo(T) could somehow become measurable in 
an appropriate experiment. One could naively imag- 
ine that it would suffice to investigate very high sound 
frequencies, as for example by going to sufficiently high 
momentum exchange Q in an x-ray Brillouin scattering 
experiment, to eliminate the viscous relaxation term in 
(2b) by forcing f2T t h 3> 1- Consider however the value 
T t h — 5 x 10~ 14 s which is found at room T in Fig. 6. 
To achieve firth — 1 would require Q/2ir ~ 3 THz which 
is already above uj c /2tt ~ 2 THz. One runs then into 
a regime in which the sound waves are presumably hy- 
bridizing j2(J with the Si04 librational mode j^Jj, pro- 
ducing the boson peak 30]. Although such a measure- 
ment might give an indication about Voo, it is certainly 
not a clean cut situation. The other alternative would 
be to reduce the temperature so as to increase the value 
of r t h ■ This also soon meets its limits as the size of the 
viscous relaxation becomes negligible below ~ 100 K, as 
seen from Fig. 7. We thus find that it would unfortu- 
nately be difficult to obtain a direct measurement of the 
unrelaxed sound velocity. 

To conclude, it seems that we have for the first time 
found a glass in which the Brillouin sound damping ap- 
pears to be entirely controlled by the network viscosity. 
This already suggests that it is probably not appropriate 
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to simply n egle ct the anharmonic damping as proposed 
for silica in [33. This illustrates that the question is of 
importance, as it finally ties in with much more funda- 
mental and debated issues, like the origin of the thermal 
anomalies of glasses 4]. We are fully aware that the noise 
in Fig. 2b is still far from satisfactory. We could only 
progress on this issue if we had an appreciably more ho- 



mogeneous sample than the already excellent one used 
in this study. If so, it would be relatively straightfor- 
ward to increase the spectrometer performance by using 
a spherical Fabry-Perot of larger spacing, such as 10 or 
25 cm, the latter giving a factor five improvement over 
the present resolution. 
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